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Abstract 

Q_i! In color dipole gBFKL dynamics, we describe the emerging gBFKL phe- 

nomenology of a subasymptotic energy dependence of the diffraction slope 
and discuss possibilities of testing the gBFKL predictions in exclusive photo- 
and electroproduction of vector mesons V at HERA. A substantial shrinkage 
of the diffraction cone j*p — > Vp processes from the CERN/FNAL to HERA 
range of energy W is predicted. This subasymptotic shrinkage is faster than 
expected from the small slope of pomeron's Regge trajectory a'jp. We point 
out that the diffraction slope is a scaling function of (my + Q 2 ), what relates 
production of different vector mesons. 
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Is the QCD pomeron a fixed or moving singularity? Can the two options be distinguished 
experimentally in hard diffraction processes at HERA? These pressing issues are adressed 
in this paper to. 

The early discussion on the BFKL (Balitsky-Fadin-Kuraev-Lipatov [1]) pomeron cen- 
tered on the scaling a$ = const approximation with infinite propagation radius R c of 
(massless) gluons. In this approximation, for a lack of the dimensional scale, the BFKL 
pomeron is a fixed cut in the j-plane: — oo < j < ajp(t) = ajp(0) = 1 + An>, which also 
implies the energy-independent diffraction slope at t — 0. The scaling approximation is 
not self-contained, though, because the diffusion BFKL Green's function makes high-energy 
behavior even at short distances increasingly sensitive to the nonperturbative large- distance 
contribution [1-3]. The introduction of running QCD coupling as{r) into the generalized 
BFKL (gBFKL) equation for color dipole cross section [3] further amplifies an intrusion of 
large-distance effects. Furthermore, in our Ref. [4] we have shown that breaking of scale 
invariance by a running as(r) coplemented by the finite gluon propagation radius R c [3], 
profoundly changes the very nature of the gBFKL pomeron from a fixed cut of the scaling 
as = const approximation to a moving cut with the finite Regge slope a'jp of the pomeron 
trajectory. This conclusion [4] derives from the fact that gBFKL diffraction slope B(£,r) 
for scattering of a color dipole of size r has the asymptotic Regge growth 

S(e,r) = S(0,r) + 2aU, (1) 

where £ = \og(W 2 /s Q ), W is the c.m.s. energy and s = m 2 in hadronic scattering, s = Q 2 
and £ = log(-) in deep inelastic scattering and in leptoproduction of vector mesons so = 
my + Q 2 . In the gBFKL dynamics [3-5] a dimensionful a' w is a noneperturbative quantity 
related to the nonperturbative infrared parameter of the model - the gluon propagation 
radius R c . For the preferred R c = 0.27 fm, quite a small a'jp = 0.072 GeV~ 2 is found 
(the Regge phenomenology of soft hadronic scattering suggests [] ot' so ^ t ~ 0.2GeV~ 2 [6,7]). 
A strong argument [3-5,8-10] for using running a>s(r) is that only in this case the short- 
distance (Leading-LogQ 2 ) limit of the gBFKL equation matches the GLDAP equation [11]. 

^^Here one must bear in mind that because of the so-called absorption corrections the observed &' so f t is 
typically larger that the Regge slope a' e ^ in the bare pomeron amplitude [7] . 
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A dramatic impact of scale invariance breaking on the intercept Ajp found in [3] (for earlier 
work see [9]) has recently been confirmed in [10,12] using a very different technique. 

How the finding that the gBFKL pomeron is a moving cut can be tested experimentally? 
Based on a solution of the gBFKL equation [4] for B(£,r), in this paper we discuss the 
emerging gBFKL phenomenology of the diffraction slope. First, with plausible and partly 
tested boundary conditions, we show that at attainable subasymptotic energies (comprising 
the HERA energies), the local Regge slope a' e jf(£, r) = \dB(£, r) / <9£ is much larger than the 
avove cited a' w and is surprisingly close to Oi' so f t ~ 0.2 GeV -2 . Second, we find £?(£, r) of a 
magnitude which is intriguingly close to the experimental determinations. Third, we discuss 
scaling properties of B(£,r) which relate diffraction slope for different exclusive reactions 
7*p — ► Vp (V = p°,uj°, 0°, J/^, T, ..) and can be tested at HERA. Here we focus on the 
fact [13,14] that exclusive vector meson production probes the dipole cross section <r(£,r), 
and B(£,r) and a' e ff(£,r) thereof, at r ~ rs, where the scanning radius 

r s « 6 . (2) 
^Q 2 + K 

We recall first the gBFKL equation [3,5] for color dipole cross section <r(£,r): 



= /C®<7(r,0 = 



^ f j2 -> 2 



.2 

gs(Ri)Ki(PGfh)— - gs{R2)K l {fi G p2)— 

Pi P2 



[ofori + vfori-vfor)] . (3) 
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Here the kernel K, is related to the flux of Weizsacker- Williams (WW) soft gluons generated 
by the q-q color dipole source, r is the q-q separation and pi j2 are the q-g and q-g separations 
in the two-dimensional impact parameter plane, K v (x) is the modified Bessel function, 
£{p) = /iGgs(p)Ki(p>GP)^ = —Qs(p)V pK (ficp) describes a Yukawa screened transverse 
chromoelectric field of the relativistic quark and 



Pg 



9 s{Ri) K\{pgPi)— - gs(R2)K 1 ( P L G p 2 )^ 



2 

2 



= \Efc)-E(fh)\ 2 (4) 



Pi P2 

gives a flux (a modulus of the Poynting vector) of WW soft gluons in the qqg state and 
§[<t(£, Pi)+<j(£, p2)—<j{i, r)\ is a change of cross section for the presence of the WW gluon [5]. 



The running QCD charge gs(f) = \/^ a s( r ) must be taken at the shortest relevant distance 



Ri = min{r,pi} and tn the numerical analysis [3] an infrared freezing a,s{q 2 ) < 0.82 has 
been imposed on the three-flavor, one-loop as(q 2 ) = 47r/[9 log(g 2 /A 2 )] with A = 0.3 GeV. 
The preferred choice R c = 0.27 fm gives Ajp = 0.4 [4] and leads to a very good description 
[15] of the HERA data on the proton structure function at small x. 

The nonperturbative infrared parameter R c = 1/pc has a very clear meaning of a cor- 
relation (propagation, Ukawa screening) radius for perturbative gluons. The above infrared 
regularization is not unique, but the crucial transversality property of WW gluons holds 
independent of the specific model for screening and only the numerical results can some- 
what change. At r, pi )2 <C R c and in the as = const approximation, the scaling BFKL 
equation is obtained ([3,5], see also Mueller and Patel [16]), in the Leading-LogQ 2 regime 
of r < p 1>2 < R c , the gBFKL Eq. (g) matches [3,5] the GLDAP equation (see also [9,10]). 

Generalization of (|[) to the equation for diffraction slope B(^,r) proceeds as follows [4]: 
If the impact-parameter representation, cr(£,r) = 2J d 2 bT(C,,r,b) and the diffraction slope 
B(C,r) at t = equals B(£,r) = |(6 2 ) = A(£,r)/<r(£,r), where 

A(£,r)= Jd 2 bb 2 T(^r,b) : (5) 

— # — # 

r(£, r, b) is the profile function and b is the impact parameter defined with respect to the 
center of the q-q dipole. In the qqg state, the qg and qg dipoles have the impact parameter 
b + \p%\ and, undoing the impact parameter integrations in (§), one finds: 



/C®r(£,r,6) = — J d 2 p 1 p 2 G 



g s {Ri)Ki(p G pi)— - g s (R 2 )K x {p jG p 2 ) — 
Pi P2 



2 



x [r(e, pi, b + l -p 2 ) + r(£, P2 , b + i#) - r(e, r, b)} . (6) 



It is convenient to separate from -B(£, r) a purely geometrical component |r 2 due to the color 
dipole's elastic form factor and consider //(£, r) = A(£, r) — |r 2 cr(£, r) . Then, the calculation 
of the moment (|5]) in (^) leads to an inhomogeneous equation 



dr)(£,r) 3 



J d 2 p 1 p 2 G 



9s(R\)K\{tJ>GPi)— - 9s{R2)K 1 (p G p 2 ) — 

Pi P2 



2 



d£ 8tt 3 

x [vtt, Pi) + V(C, P2) - ri(S, r) + - {pi + p\ - r 2 ) [a(p 2 , + 

= /C®tj(£, 0+0(60, (7) 
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g s (Ri)Ki(n G pi)— - gs{.R2)Ki(ii G p 2 ) — 

Pi P2 



/3(£, r) = C ® r) = ^sf d 2 p, p 2 G 

x(^ + P^-r 2 )[ ( x(p 2 ,0 + ^(pi,0] • (8) 



A homogeneous Eq. ([7]) coincides with Eq. ([$]), which has several consequences. First, if 
0"i (£, r ) and 771 (£, r) are solutions of Eqs. (|) , (0) , then 772 (£,?") = rji (£, r) + A6 ■ 0"i (£, r) with 
Ab = const also is a solution of Eq. (^) with £? 2 (£,?") = r) + Ab. Second, the Regge 
growth of B(£,r) can be driven only by the inhomegeneous term /3(£,r). 

The detailed arguments for a'jp 7^ 0, based on properties of eigenvalues of the gBFKL 
equation (Q) and of the kernel C have been presented in [4] and need not be repeated here. 
We only cite the order of magnitude estimate 

3 f ,n _ / \ 9 9 „Q / •. 3 



lo7r z 



■ J d 2 fa s {r)p 2 G r 2 Kl{p G r) ~ — i? c 2 a 5 (/? c ) , (9) 



which clearly shows the connection between the dimensionful a'jp and the nonperturbative 
infrared parameter R c . We determine a^p from the numerical solution of Eqs. as the 

£ — > 00 limit of «gjj(£,r). The resulting i? c dependence of a'jp is shown in Fig. 1 and is 
very steep in the opposite to a weak independence of the intercept Ajp found in [3] . 

The gBFKL dipole cross section cr(£, r) sums the Leading-Log(W /2 ) multigluon produc- 
tion cross sections. Consequently, as a realistic boundary condition for the gBFKL dynamics 
one can take the two-gluon exchange Born amplitude 

16tt f 2r a s (k 2 )a s (K 2 ) 



l07T f -* 

W(£ ,r,g) = — J d 2 k 



[{k + W + vMh-W + pl] 
x [cos(igr) - cos(fcf)] • [G^qj - G 2 {k + ~£ -£ + -£)] . (10) 

Here we use the normalization IuL4(£ ,r, g = 0) = er(£o,r), k, 2 = min{A; 2 , C 2 r~ 2 }), C ~ 1.5 
[17], k ± |<f are the momenta of exchanged gluons, G\(q) and (^2(^1,^2) are the single- 
and two-quark form factors of the proton probed by gluons. The former is customarily 
approximated by, and the latter in simple quark models can be related to, the charge form 
factor of the proton (for instance, see [18,19]). A very good quantitative gBFKL description 
of the HERA data on the small-x proton structure functions is obtained [15] if the Born 
cross section Eq. ( |10"D with R c = 0.27 fm is taken as a boundary condition for the gBFKL 
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equation (|3|) at the Bjorken variable xo = 0.03, i.e., at £ = log — = 3.5. Hereafter we only 
consider R c = 0.27 fm. 

The Born approximation for B(£ ,r), shown in Fig. 2, has nice properties which admit 
a simple interpretation. First, for r < R c , the Born amplitude (|1(J) is dominated by pertur- 
bative k 2 ~ r~ 2 > jj? G for which G2(k, —k) vanishes. Consequently, A(^,r,q) oc G\(q) and 
B(£ ,r <C R c ) = |(_R 2 ) p , which is a generic perturbative result and holds beyond the above 
simplified derivation. Second, the r dependence of the Born approximation B(C,q, t) can be 
cast in a very symmetric and intuitively appealing form 

B(C , r) = ±(Rj) + l -r 2 + 5B = l -(r 2 cms ) beam + \{r 2 cms ) target + 5B . (11) 

The found small departure from the dipole size dominated slope (|TTD, \SB\ < 1 GeV~ 2 , 
is evidently related to the small R 2 ~ 2 GeV~ 2 <C ^(R-lh)- O ne can ar gue that for the 
same reason the r dependence of B(^o,r) must be insensitive to details of the infrared 
regularization. The negligible 5B in Eq. ( |TTD implies a simple boundary condition 77 (£o, r ) ~ 
-B(£o? O)cr(^O) r ) and an energy independent contribution ps _B(£ ,0) to the diffraction slope 
£>(£,r). Consequently, the energy dependence of B(£,r) and a' e ^(C,,r) are governed by the 
inhomogeneous term /3(£,r) of Eq. (H) and only depend on the input dipole cross section 
(|10|) , which already has succesfully been tested [15] against the DIS data. The gluon- 
probed radius of the proton R± is a phenomenological parameter to be determined from the 
experiment. For the evaluation purposes, we approximate R\ by the proton charge radius 

Rch'- 

fl(£o, r « R c ) « l -{R 2 ch ) p « 5.8 GeV- 2 . (12) 

In Fig. 3 we present the effective Regge slope a' e jj(^,r) which follows from the solution 
of coupled Eqs. (0,f?U§) with the above described boundary conditions. At £ — > 00, a' e jj-(£,,r) 
tends to a r-independent a'jp = 0.072 GeV -2 . The very slow onset of the limiting value 
a'jp = 0.072 GeV -2 correlates nicely with the very slow onset of the gBFKL asymptotics of 
cr(£,r) and of the proton structure function [3,8,15]. Slight oscillations in a^(£, r) at small 
£ are a consequence of the oscillatory r-dependence of large-i/ eigenf unctions E(v, r) (for 
the related discussion see [3]). 
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The above results for a' e jj(l;, r) can be tested in exclusive diffractive processes 7*p — > Vp. 
At fixed £, the dependence on the vector meson mass my and the photon virtuality Q 2 is 
concentrated in Eq. (|2]) for the scanning radius r$ and we predict the same B(£,rs) and 
a e//(£; r$) for all reactions with the same £ and rs and/or {m v + Q 2 ), provided that rs < Rv 
[14,20]. For instance, r 5 (T, Q 2 = 0) « r s (J/tf , Q 2 ~ 120 GeV 2 ) « r 5 (p°, Q 2 ~ 200 GeV 2 ) « 
0.13 fm; r 5 (J/^,Q 2 ~ 30 GeV 2 ) « r 5 (p°,Q 2 ~ 60 GeV 2 ) « 0.2 fm; r 5 (J/^,Q 2 = 0) « 
r S (p°,Q 2 ~ 20 GeV 2 ) w 0.4fm and r 5 (p°,Q 2 ~ 3.5 GeV 2 ) « 0.76 fm (we refer to the 
dominant production of the longitudinally polarized p° [14]). The hatched areas in Fig. 3 
indicate a variation of a' e jj(C,,r) over the HERA range of c.m.s energy W = 50 — 200 GeV 
(the higher values of a e ff(£,r) correspond to a lower W). Remarkably, at the HERA and 
lower energies, the subasymptotic a e ff(£,r) ~(0. 15-0.2) GeV -2 is quite large, close to ot' so ^ t 
and by the factor 2-3 larger than the asymptotic value a' w = 0.072 GeV -2 . 

The results for the diffraction slope in an approximation (|P2|) are presented in Fig. 4. 
According to Eqs. (PiPD, the diffraction slope decreases with Q 2 and the vector meson mass. 
The absolute value of B(£,r) depends on the assumed gluon probed radius of the proton 
R±, whereas the rate of the shrinkage a' e ^{^,r) does not. The experimental determination 
of [R\) is of great interest by its own. For all the exclusive 7*p — > Vp reactions we find a 
substantial rise of B(£,rg), by ~ 1.5 GeV -2 , from the CERN/FNAL fixed target range of 
W ~ 10 GeV to the HERA collider energy W ~ 100 GeV. Here we have taken R c = 0.27 fm, 
the experimental measurement of this rise would be the best constraint on R c . 

The above analysis refers to the (infrared regularized) perturbative gBFKL scattering 
amplitude. For the description of semi-perturbative and soft scattering phenomena, one 
needs to complement the gBFKL dipole cross section by the soft nonperturbative, non- 
gBFKL, component cr^ npt \^, r), which at r > R c overwhelms the above discussed perturba- 
tive gBFKL cross secrion cr^(£,r) [14,15]. This non-gBFKL soft cross section has only a 
marginal effect on the proton structure function at large Q 2 [15], but contributes substan- 
tially to vector meson production amplitudes unless rs < R c [14]. Of the above 7*p — * Vp 
reactions, only a real (and a virtual) photoproduction of the T and a virtual photoproduc- 
tion of the J/W at Q 2 > 20 GeV 2 and of the p° at Q 2 >40GeV 2 are purely perturbative 
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processes. For instance, in real photoproduction of the J/^f, the nonperturbative soft con- 
tribution makes ~ 50% of the photoproduction amplitude. Although the mechanism of this 
nonperturbative soft interaction is not well undesrtood, the geometrical size dependence 
( |TTD could roughly be applicable also to a soft production mechanism (for an example of 
the nonperturbative model see [21]). Furthermore, because we find a' e ^ ~ a ' so fv as a P oor 
man's approximation, we can use the above gBFKL evaluations of the slope for the soft 
scattering too. The estimates of B(£,rs) for a semi-perturbative real photoproduction of 
the J/ty and moderately virtual photoproduction of the p°, shown in Fig. 4, assume this 
poor man's approximation. 

The available data on the J/ty photoproduction are not yet conclusive. A direct selection 
of purely elastic events was performed only in the FNAL E401 experiment [22] with the result 
B(J/^;Q 2 = 0) = 5.6 ± 1.2 GeV -2 . The CERRN NMC result, B(J/^;Q 2 = 0) = 5.0 ± 
1.1 GeV -2 , comes from a sample which contains background excitation of the target proton 
[23]. These values of B(J/^;Q 2 = 0) for W ~ 15 GeV are consistent with our estimates 
shown in Fig. 4. At the present stage of HERA experiments, a direct rejection of proton 
excitation is not yet possible [24,25]. The model-dependent analysis of the first HI data at 
(W) = 90 GeV gives estimates of B(J/^;Q 2 = 0) ranging from 8.1 to 4.9 GeV" 2 [24], the 
first ZEUS data give B{J/^\ Q 2 = 0) = 4.5 ±1.4 GeV" 2 [25]. One needs the higher accuracy 
data from the both fixed target FNAL/CERN and the collider HERA experiments. Here 
one must bear in mind a well known rapid rise of the diffraction slope towards t = [26]. We 
calcuate the diffraction slope at t — 0, whereas the experimental data on the vector meson 
production correspond to a slope evaluated over quite a broad range of t, typically up to \t\ < 
1 GeV" 2 , which may underestimate B(V; Q 2 ). The ZEUS data on virtual photoproduction 
of the p° mesons give B(p°; 7 < Q 2 < 25 GeV 2 ) = 5.1 + 1.2 - 0.9(stat) ± l.O(syst) GeV" 2 
[27], which is close to the above cited B(J/^/; Q 2 = 0), in agreement with our (rriy + Q 2 ) 
scaling. This result is also consistent with the predicted decrease of B(p°; Q 2 ) from real 
photoproduction in which B(p°; Q 2 = 0) = 10.6 ± OA(stat) ± 1.5(syst) GeV" 2 [28]. 

To summarize, the purpose of the present paper has been an exploration of the phe- 
nomenology of forward diffractive scattering which emerges from the gBFKL dynamics. 
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The intrusion of large distance effects leads to a shrinkage of the diffraction cone driven by 
the nonperturbative slope of the pomeron's Regge trajectory a'jp. We presented the first 
evaluation of the energy dependence of the gBFKL diffraction slope. An interesting finding 
is a large subasymptotic value of the effective Regge slope a' e jf(£,r), which is by the factor 
~(2-3) larger than a'jp. For exclusive production of all the vector mesons at HERA, we find 
a substantial, by ~ f .5 GeV~ 2 , rise of the diffraction slope -8(7* — > V) from from the fixed 
target CERN/FNAL to the collider HERA energy. From the scanning property in vector 
meson production, we predict the scaling dependence on the variable (my + Q 2 ), which 
allows to relate diffracton slope for different reactions. Strong dependence of a'jp on the 
gluon propagation radius R c makes the rise of the diffraction slope a sensitive probe of R c . 
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Figure captions: 

Fig. 1 - The Regge slope a'jp of the pomeron trajectory vs. the propagation radius of 
gluons R c . 

Fig. 2 - The solid line shows the diffraction slope B(^ Q , r) for color dipole-proton scat- 
tering in the two-gluon exchange Born approximation. The dashed straight line shows the 
geometrical law Eq. ( |TTD with. The difference between the dashed line and solid curve de- 
scribes 5B. 

Fig. 3 - The energy (£) dependence of the effective Regge slope a' e ^{^,r) vs. £ = 
log(l/x e //) = \og\W 2 /{Q 2 + my)] changing by A£ = 2 from £ = 5.5 (the curve (a)) 
to £ = 13.5 (the curve (e)). The curve (/) is for £ = 17.5. For reference, the curves 
(a), (b), (c), (d), (e) and (/) correspond to x eff = W 2 /(Q 2 + m 2 v ) = 5 ■ 1(T 3 , 7 • 10" 4 , 9 ■ 
10 -5 , 1.25 • 10~ 5 ,1.6 • 10~ 7 and 3 • 10~ 8 , respectively. The horizontal line (g) shows the 
asymptotic value a' w = 0.072 GeV~ 2 . The hatched areas A — D show the range of vari- 
ation of a' e ^{^,r s ) over the HERA energy range W =(50-200) GeV at a scanning radius 
r$ relevant to the reaction 7*p — ► Vp at the photon's virtuality Q 2 : (A) rs = 0.12fm, 
T(Q 2 ~ 0), (B) r s = 0.21 fm, J/^(Q 2 ~ 30 GeV 2 ), (C) r s = 0.4 fm, J/V(Q 2 ~ 0) and 
p°{Q 2 ~ 20 GeV 2 ), (D) r s = 0.76 fm, p°{Q 2 ~ 3.5 GeV 2 ). 

Fig. 4 - The c.m.s. energy W dependence of the diffraction slope B(£, rs) at a dipole size 
(scanning radius) rg relevant to different diffractive ^*p — > Vp processes: (a) rs = 0.12fm, 
T(Q 2 ~ 0), (b) r s = 0.21 fin, J/^(Q 2 ~ 30 GeV 2 ), (c) r s = 0.4 fm, J/^(Q 2 ~ 0) and 
p°(Q 2 ~ 20 GeV 2 ), (d) r 5 = 0.76 fm, p°{Q 2 ~ 3.5 GeV 2 ). 
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